Universality behind Basquin's law of fatigue 
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One of the most important scaling laws of time dependent fracture is Basquin's law of fatigue, 
namely, that the lifetime of the system increases as a power law with decreasing external load 
amplitude, tf ~ o"o" a , where the exponent a has a strong material dependence. We show that in 
spite of the broad scatter of the Basquin exponent a, the fatigue fracture of heterogeneous materials 
exhibits intriguing universal features. Based on stochastic fracture models we propose a generic 
scaling form for the macroscopic deformation and show that at the fatigue limit the system undergoes 
a continuous phase transition when changing the external load. On the microlevel, the fatigue 
fracture proceeds in bursts characterized by universal power law distributions. We demonstrate 
that in a range of systems, including deformation of asphalt, a realistic model of deformation, and 
a fiber bundle model, the system dependent details are contained in Basquin's exponent for time to 
failure, and once this is taken into account, remaining features of failure are universal. 

PACS numbers: 46.50.+a, 62.20.Mk, 61.82.Pv 
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Disordered media subject to sub-critical external loads 
present a time dependent macroscopic response and typ- 
ically fail after a finite time Such time dependent 
fracture evidently plays a crucial role in a large variety 
of physical, biological, and geological systems, such as the 
rupture of adhesion clusters of cells in biomaterials un- 
der external stimuli Q , the sub-critical crackgrowth due 
to thermal activation of crack nucleation 0, Q , creep Q 
and fatigue fracture of materials @,0|j and the emergence 
of earthquake sequences [f| . One of the most important 
scaling laws of time dependent fracture is the empirical 
Basquin law of fatigue which states that the lifetime tf of 
samples increases as a power law when the external load 
amplitude cr decreases, tf ~ a^ a Q. The measured 
values of the Basquin exponent a typically vary over a 
broad range indicating a strong dependence on material 
properties [H, [l(| 11 1 . 



In this Letter we study the fatigue fracture of hetero- 
geneous materials focusing on the underlying microscopic 
mechanism of the fatigue process and its relation to the 
macroscopic time evolution. We develop two generic 
models of time dependent fracture, namely, a fiber bundle 
model and a discrete element approach, which both cap- 
ture the most important ingredients of the fatigue failure 
of disordered materials. Analytic solutions and computer 
simulations reveal that the models recover the Basquin 
law of fatigue, whose exponent is determined by the dam- 
age process. We show that, as a consequence of healing, 
a finite fatigue limit emerges at which the system under- 
goes a continuous phase transition from a regime where 



macroscopic failure occurs at a finite time to another one 
exhibiting only partial failure in the system having an in- 
finite lifetime. Based on analytic solutions, we propose a 
generic scaling form for the macroscopic deformation. On 
the microlevel the fatigue of the material is accompanied 
by an avalanche activity where bursts of local breakings 
arc triggered by damage sequences. Wc demonstrate an- 
alytically that the microscopic bursting activity underly- 
ing fatigue fracture is characterized by universal power 
law distributions which implies that the non-universality 
of the Basquin exponent at the macro-level is solely due 
to the specific degradation process of the material. 

First we consider a mean field model of fatigue frac- 
ture, namely, a fiber bundle model (FBM) where fibers 
fail either due to immediate breaking or to ageing 14 1. 



For the load redistribution after failure events, equal load 
sharing is assumed so that all the fibers carry the same 
load [12(. During the evolution of the system, a fiber 
breaks instantaneously at time t when the load on it p(t) 
exceeds the local tensile strength p\ h (i = 1, . . . , N). All 
intact fibers accumulate damage c(t) due to the load p(t) 
that they have experienced and break when c(t) exceeds 
the local damage threshold c\ h (i = 1, . . . , N). The ac- 
cumulated damage c(t) up to time t is obtained by in- 
tegrating over the entire loading history of the specimen 



c(t) = afe' 



~p(t')' y dt', where a > is a scale pa- 



rameter, while the exponent 7 > controls the rate of 
damage accumulation [h], 11 1. To capture damage re- 
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covery in the model due to healing of microcracks 
thermally activated rebinding of failed contacts @, 
introduce a memory term in the above damage law of ex- 
ponential form whose characteristic time scale r defines 
the memory range of the system 0, 0, Q • Hence, during 
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FIG. 1: (Color online) Scaling plot of deformation-time 
curves measured experimentally on asphalt specimens [l4| . 
The scaling function S, obtained by rescaling the two axis, 
has a power law dependence on the time-to-failure. 



the time evolution of the bundle, the damage accumu- 
lated over the time interval t' < (t — r) heals. Assuming 
independence of the two breaking thresholds p t h and Cth, 
the macroscopic evolution of the system under a constant 
external load <tq can be cast into the form 



a = [l-F(c(t))] [1-G(p(t))]p(t), 



(1) 



where G and F denote the cumulative distributions of 
Pth and c t h, respectively. We solved Eq. (Q} analytically 
obtaining the load p(t) on the intact fibers at a constant 
external load Co < cr c , with the initial condition p(t = 
0) = POi where po denotes the solution of the constitutive 
equation <tq = [1 — G(po)]po [12] • Here a c denotes the 
ultimate strength of the material. The most important 
input parameters of the model calculations are a, 7 and 
t, which govern the damage accumulation. 

On the macrolevel the process of fatigue is character- 
ized by the evolution of deformation e(t) of the specimen, 
which is related to p(t) as p(t) = Ee(t), where E = 1 
is the Young modulus of fibers. Neglecting immediate 
breaking and healing, Eq. |(TJ) can be transformed into a 
differential equation for the number Nj, of broken fibers 
as dN b /dt = a/(c(t))pTJV, where /(c) = dF/dc. Us- 
ing p(t) = N<j / (N — Nb(t)), for uniformly distributed 
threshold values the exact solution of the equation of mo- 
tion Eq. (JTJ) reads 



£ (t) = <7 [(i / -i)A/]" 1/(1+7) and*/ 



o(l + 7) 



(2) 



where t / denotes the lifetime of the system. Equation pj) 
shows that damage accumulation leads to a finite time 
singularity where the deformation e(t) of the system has 
a power law divergence with an exponent determined by 
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FIG. 2: (Color online) Characteristic time scales t r and tj 
of the system. The complete FBM corresponding to Eq. ([T]) 
which includes immediate breaking and healing is solved nu- 
merically. For <r > (J; we see Basquin's law and both models 
provide a very good fit of the lifetime data of asphalt. The 
fatigue limit oi is indicated by the vertical dashed line. 



7. It is important to emphasize that tf has a power law 
dependence on the external load in agreement with 
Basquin's law of fat igue fo und experimentally in a broad 
class of materials 0, (TT 
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The Basquin exponent of the 
model therefore coincides with that of the microscopic 
degradation law a = 7. Another interesting outcome of 
the derivation is that the macroscopic deformation e(t) of 
a specimen undergoing fatigue fracture obeys the generic 
scaling form e(t) = o-qS^ctq), where the scaling function 
S has the property S[taf) - (t a - taH)~ 1/( - 1+ ^ , with 
t a = 0(1 + 7) and the scaling exponents are S = 1 and 
(3 = 7. Figure [1] presents a verification of this scaling law 
on experimental data from asphalt specimens obtained 
at two different load values [1J|. The good quality data 
collapse obtained by rescaling the two axis and the power 
law behavior of 5 as a function of the time-to-failure 
demonstrates the validity of our scaling relation. 

Healing dominates if for a fixed load ctq the memory 
time r is smaller than the lifetime obtained without heal- 
nl g r ^ tf{°~o, T = +00). Then, a threshold load ai 
emerges below which the system relaxes, i.e., the de- 
formation e(t) converges to a limit value with a char- 
acteristic relaxation time t r resulting in an infinite life- 
time. Figure [2] presents the characteristic time scale of 
the system varying the external load over a broad range. 
The results from numerical simulations with the com- 
plete FBM (i.e., including immediate breaking and heal- 
ing) are in excellent agreement with the measured life- 
time of asphalt samples for er > recovering also the 
Basquin exponent [14| . The regime below ai is of par- 
ticular importance in gcodynamics where memory effects 
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FIG. 3: The relaxation time t r and lifetime tf as a function 
of \ai — CTo| for different disorder distributions (uniform and 
Weibull) and 7 exponents. The straight lines have slopes 
-1/3 and -2/3. 

take place during cyclic loading of rocks with a stress am- 
plitude increasing from one cycle to the next [lq . It is 
important to note that approaching the fatigue limit a\ 
from either side, the characteristic time scale diverges. 
Figure [3] shows that both the relaxation time t r and the 
lifetime tf follow a power law as a function of the dif- 
ference from the fatigue limit with distinct exponents: 
t r ~ (<?7 — oo) -1 ^ 3 and tf ~ (er — &i)~ 2 ^- We stress that 
the exponents neither depend on the disorder distribu- 
tions (F and G) nor on the details of the damage law (a, 
7 and t), i.e., they are universal implying a continuous 
phase transition at the fatigue limit 07 between partial 
failure and macroscopic fracture (see Fig. [3]) . 

Our calculations revealed that the Basquin law of life- 
time emerges on the macrolevel as a consequence of the 
competition between the two microscopic failure mecha- 
nisms of fibers. Rewriting Eq. (TTJ) in the form of the con- 
stitutive equation of simple FBMs as <tq/ [1 — F(c(t))] = 
[1 — G(p(t))] p(t) it can be seen that the slow damage 
process on the left hand side quasi-statically increases the 
load on the system: ageing fibers accumulate damage and 
break slowly one-by-one in the increasing order of their 
damage thresholds c\ h . After a number of damage 
breakings, the emerging load increment on the remaining 
intact fibers can trigger a burst of immediate breakings. 
Since load redistribution and immediate breaking occur 
on a much shorter time scale than damage accumulation, 
the entire fatigue process can be viewed on the microlevel 
as a sequence of bursts of immediate breakings triggered 
by a series of damage events happening during waiting 
times T, i.e., the time intervals between the bursts. The 
microscopic failure process is characterized by the size 
distribution of bursts P(A), damage sequences P(Ad), 



and by the distribution of waiting times P(T). At small 
loads (To cr c most of the fibers break in long damage 
sequences, because the resulting load increments do not 
suffice to trigger bursts. Consequently, the burst size dis- 
tribution P(A) has a rapid exponential decay. Increas- 
ing Co the total number of bursts rib increases linearly 
rib ~ fo and a power law regime of burst sizes emerges 
P(A) ~ A~^ with the well-known mean field exponent 
of FBM £ = 5/2 [l5j]. When macroscopic failure is ap- 
proached (To a c the failure process accelerates such 
that the size Ad and duration T of damage sequences de- 
crease, while they trigger bursts of larger sizes A, and fi- 
nally macroscopic failure occurs as a catastrophic burst of 
immediate failures. Since in the limiting case of Co — * °c 
a large number of weak fibers breaks in the initial burst, 
we found that the distribution -P(A) has a crossover to 
a smaller exponent £ = 3/2, in agreement with Ref. [l5| . 
After the linear increase, the number of bursts rib has a 
maximum at ao/a c ~ 0.4 and rapidly decreases to 1 as 
<7 C is approached. All these results are independent of 7, 
a, and r. 

Since damage events increase the load on the remain- 
ing intact fibers until an immediate breaking is triggered, 
the size of damage sequences A^ is independent of the 
damage characteristics c(t) and F(cth) of the material, 
instead, it is determined by the load bearing strength dis- 
tribution G(j>th) of fibers. Under broad conditions this 
mechanism leads to an universal power law form with 
an exponential cutoff P(A<j) ~ A7 1 exp(— Ad/ (Ad)), 
where (A^) ~ cr^ 1 . The damage law c(t) of the ma- 
terial controls the time scale of the process of fatigue 
fracture through the temporal sequence of single damage 
events. In damage sequences fibers break in the increas- 
ing order of their damage thresholds c\ h which determine 
the time intervals At between consecutive fiber break- 
ings. Analytic calculations showed that P(At) has an 
explicit dependence on 7 as P(At) ~ At - ' 1-1 / 7 ), how- 
ever, the duration of sequences T = J2f=i ^j, i.e., the 
waiting times between bursts follow an universal power 
law distribution P(T) ~ T^cxp^T/ (T», where only 
the cutoff has 7-dependence (T) ~ c ^ 1+7 ' 1 (see Fig. SJ. 

The macroscopic lifetime tf of a finite system can be 
related to characteristic quantities of the microscopic fail- 
ure process as tf — Yl7=i Ti, from which the average life- 
time can be obtained in the form tf « (rib) (T). In the 
load regime where the generic scaling laws of the distri- 
butions P(A), P(A d ), and P(T) prevail, this leads to 
the form tf ~ er 7 in agreement with the Basquin law 
Eq. ^) of the system. The results demonstrate that the 
Basquin law of lifetime on the macro-scale is a fingerprint 
of the scale-free microscopic bursting activity, with the 
material dependence entering only through the damage 
law determining the waiting times between bursts. Ex- 
perimentally, the microscopic fracture process underlying 
fatigue can be monitored by the acoustic emission tech- 
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FIG. 4: (Color online) Rescaled plot of waiting time distri- 
butions P(T) obtained by FBM simulations. Inset: Corre- 
sponding results of two-dimensional DEM simulations. 



nique and by direct optical observations [j], 0, [l^, 11 1. 
Sub-critical cracking has recently been found to produce 
a power law distribution of step sizes of the advancing 
crack in agreement with our predictions on the size dis- 
tribution of bursts [H . 

In order to study the effect of stress concentration 
and crack growth in fatigue fracture, we also developed 
a discrete element model (DEM) [13[ in which we dis- 
cretize a two-dimensional disc-shaped specimen in terms 
of randomly shaped convex polygons connected by elas- 
tic beams. The beams fail either due to immediate 
breaking or damage which are coupled in a single fail- 
ure variable q{t) = p(t) + a f* e -(*-*')/ r p(t') 7 ^'- Here 
p{t) describes the deformation state of the beam tak- 
ing into account both stretching and bending p(t) = 
(e/e t h) 2 +max(|9i|, |@2 being e the longitudinal 
deformation, 0i and O2 the bending angles at the two 
ends of the beam, and e t h and <d t h denote the threshold 
values a beam can sustain under stretching and bend- 
ing, respectively. As a consequence, the parameters a, 
7, and r play the same role as their counterparts in our 
FBM. The time evolution of the system is followed by 
numerically solving the equations of motion of polygons. 
The breaking criterion q(t) > 1 is evaluated at each time 
step and beams which fulfil the condition are removed 
jl3j . We study the fatigue fracture under diametric com- 
pression of discs with constant stress <7o (Brazil test). 
Figure [5] shows that DEM provides also an excellent fit 
of the lifetime data of asphalt specimens [lij ■ DEM sim- 
ulations revealed that in the presence of stress concentra- 
tions bursts are spatially correlated and they can be iden- 
tified as sudden advancements of slowly growing cracks. 
DEM results on burst characteristics also show power law 
behavior as the mean field FBM, but with different expo- 
nents due to the two-dimensionality of the model. The lo- 



calized stress concentration built up around cracks gives 
rise to higher values of the exponents of the size distribu- 
tion of bursts P(A) ~ A~ 2 - 7 , and of damage sequences 
-P(Ad) ~ AT , while for the waiting time distribution 
P(T) the DEM exponent falls very close to the mean field 
value (see Fig. [J}. The results proved to be independent 
of the value of 7. 

Although the exponent of Basquin's law depends on 
the microscopic damage accumulation, we found an as- 
tonishing spectrum of universal features hidden behind 
this originally empirical law. On one hand we discovered 
in the experimentally relevant situation of finite dam- 
age memory a continuous phase transition between par- 
tial failure and macroscopic rupture. On the microscopic 
level of individual breaking events we showed that the 
separation of time scales of the two competing failure 
mechanisms leads to a bursting activity, where we dis- 
closed several universal scaling laws in the distributions 
and determined their exponents as well in mean-field as 
in two dimensions. In summary our approach provides 
a direct connection between the microscopic mechanisms 
constituting the main ingredients of the model (i.e., im- 
mediate breaking, damage accumulation and healing of 
microcracks) and the macroscopic behavior of the fatigue 
process. The (macroscopic) exponent from Basquin's law 
coincides with the (microscopic) exponent of the degra- 
dation law, namely = 7. Following a slightly differ- 
ent pathway, our methodology is also capable to show 
explicitly the bridge between the (universal) mechanism 
related with the scale-free bursting activity at the micro- 
scale and the (non-universal) lifetime law of the material 
at the macro-scale. 

This work opens up new experimental challenges. Our 
scaling relation of the macroscopic deformation should 
be verified on various types of materials, after which it 
could help to extract the relevant information from fa- 
tigue life measurements. For instance, it would be inter- 
esting to check our theoretical predictions with fatigue 
measurements performed at very low external loads, i.e., 
for (To ~ 07. More precisely, in the infinite lifetime limit, 
(To i$ , the experimental confirmation of the power-law 
variability with load of the relaxation time should cer- 
tainly provide some considerable insight on the role of 
healing in the entire fatigue process. For similar reasons, 
it would be also interesting to verify the distinct lifetime 
behavior obtained from the model in the other limit of 
low external loads, <tq > 07. Finally, another interesting 
outcome from our study is the statistical behavior related 
with the bursting activity during the fatigue process of a 
given material. According to our analysis, both the size 
Ad of damage sequences and magnitude T of waiting 
times between bursts should obey universal power-law 
distributions that might reflect the intrinsic features of 
the typical restructuring events taking place at the micro- 
scopic level. As a possible monitoring technique, acoustic 
emission measurements could be conducted in conjunc- 
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tion with fatigue experiments to confirm our claim for 
universality behind Basquin's law. 
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